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I .  INTRODUCTION 

Interior  ballistic  calculations  usually  require  a  knowledge  of  the 
geometry  of  propellant  grains  at  each  phase  of  burning  from  ignition  to 
burnout.  Any  calculational  scheme  for  providing  the  surface  area  or 
volume  of  a  grain  at  any  instant  is  often  referred  to  as  a  "form-function". 
Investigators  have  determined  form  functions  for  a  wide  variety  of  grain 
shapes!*2,  This  report  concerns  itself  with  providing  a  general  method 
of  constructing  form- functions  for  one  of  the  moct  popular  classes  of 
grains,  namely,  the  multiperforated  eylinu  \  .  The  shape  is  that  of  a 
right  circular  cylinder  pierced  its  entire  length  by  some  number  of 
uniform  cylindrical  perforations. 

Many  approaches  to  analyzing  particular  members  of  this  class,  for 
example,  the  seven-  or  nineteen-perforated  grains,  have  been  taken  in 
tho  past.  Instead  of  examining  many  such  special  cases  here,  a  general 
method  is  presented  that  leads  to  the  construction  of  form- functions  which 
can  determine  surface  area  and  volume  for  virtually  any  cylindrical 
grain  with  any  number  (n>2)  of  arbitrarily  placed  perforations;  tho 
special  excluded  cases  are  discussed  as  well.  The  method  neither  employs 
approximations  in  calculation  nor  does  it  assume  equal  rates  of  regression 
for  tho  inner  and  outer  grain  surfaces.  Thus  this  approach  supercedes 
several  previous  efforts  and,  moreover,  provides  tho  opportunity  to 
analyze  many  novel  cylindrical  grain  designs  in  a  straightforward 
manner.  Non-cylindrical  grains  as,  for  instance,  hexagonal-surfaced 
ones,  are  barred,  although  similar  methods  could  clearly  bo  employed. 

Finally,  two  FORTRAN  subroutines  which  perform  the  required 
computations  are  provided  and  are  used  in  connection  with  a  sample 
exercise  wherein  a  form-function  for  a  nineteen-perforated  grain  with 
three  independent  webs  is  developed. 

In  what  follows,  all  angle  measurements  are  in  radians. 


II.  OUTLINE  OF  THE  METHOD 

A.  SUBDIVISION  OF  THE  GRAIN  END  AREA 

The  circular  cross-section  of  '  lie  grain  in  question  is  first 
geometrically  decomposed  by  straight  lines  into  classes  of  mutually 
congruent  "slivers",  definod  below.  At  any  point  in  burning,  the 

1,  Corner,  J. ,  "Theory  of  the  Interior  Ballistics  of  Guns",  « John  Wiley 
S  Sons,  New  York,  1950 ,  pp.  30-55. 

2,  Stats,  J,,  "Form-Functions  for  Multicomponent  Propellant  Chargee 
Including  Inhibited  Grains  and  Sliver  Burn",  MHL  Technical  Note  371, 
Materials  Research  Laboratories ,  Maribyrnong ,  Victoria,  Australia, 
September  1975, 
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individual  contribution  of  each  class  to  the  remaining  surface  and  volume 
is  noted,  the  totals  simply  being  the  sums,  over  all  classes,  of  the 
number  of  members  in  each  class  multiplied  by  that  class's  contributions. 

There  are  two  types  of  "slivers",  hereafter  referred  to  without 
quotes:  inner  and  outer.  An  inner  sliver  is  any  region  of  cross-section 
bounded  by  a  triangle  with  the  center  of  a  perforation  at  each  vertex, 
but  with  no  portion  of  any  other  perforation  contained  within  the  boundaries. 
An  outer  sliver  is  any  region  of  cross-section  bounded  by  an  arc  of  the 
grain *s  outer  surface,  two  radii  enclosing  an  angle  less  than  n,  and  a 
chord  segment  joining  the  centers  of  two  perforations.  No  portion  of 


any  other  perforation  may  be  contained  within  its  boundaries. 

For  example,  the  grain  in  Figure  1  has  been  divided  into  two  classes 
of  slivers,  one  class  of  six  inner  slivers  and  one  class  of  six  outer 
slivers. 

A  complete  analysis  of  the  general  outer  sliver  and  the  general 
inner  sliver  follows, 

B.  OUTER  SLIVERS 

In  addition  to  the  above-mentioned  assumptions  about  the  grain  and 
its  dissection  into  slivers,  we  require  that,  throughout  burning,  the 
outer  slivers  remain  within  their  original  boundaries  until  extinction. 
Initially,  an  outer  sliver  appears  as  PQA,A,  in  Figure  2  with  A,  the 
center  of  the  grain  and  A2  and  A3  the  centers  of  the  sliver's  perforations. 
In  what  follows,  we  adhere  to  the  convention  that  side  Sj  of  triangle 
a1a2a3  *s  the  9®8ment  connecting  the  centers  of  the  sliver's  perforations. 
Further,  the  restriction  that  angle  PAjQ<ir  forces  the  grain  to  have  at 
least  three  perforations  with  not  altogether  arbitrary  placement.  We 
assume  lengths  and  S3  to  be  known. 

* 

From  the  law  of  cosines  we  have  the  following  angles 


As  burning  progresses,  the  sliver  erodes  as  shown  in  Figure  3  with, 
possibly,  some  degenerate  angles.  We  require  knowledge  of  the  grain's 
instantaneous  radius,  RAD,  grain  length,  GRL,  and  perforation  diameter, 
PRFD.  The  surfaces  of  all  perforations  are  assumed  to  erode  at  the  same 
rate,  but  this  need  not  be  the  same  rate  of  regression  as  that  of  at  the 
outer  surface.  Clearly,  burnout  of  the  sliver  occurs  when  either 
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*1 


Figure  3.  General  Outer  Sliver  During  Burning 


13 


T1  +  T2  ^  VlA2 


or 


GRL  <  0 


An  error  condition  occurs  whenever  either 


T  3  < 


A2A3A1 


or 


r4  <  A3A2A1 


at  some  step  during  burning  and  signals  violation  of  our  requirement  that 
the  sliver  stay  within  its  original  boundaries.  A  grain  exhibiting  this 
error  in  one  of  its  outer  slivers  cannot  be  analyzed  beyond  that  point 
by  the  present  method.  Very  thin  outer  slivers  are  likely  to  give  rise 
to  this  error. 

Using  the  familiar  minimum  and  maximum  functions  we  derive  by  law 
of  cosines 


T1 


T2 


T 


3 


arccos 


arc cos 


arccos 


S2"!+RAD2-l/4-PRrD2 
2'.  S'".  RAD 

S32+RAD2-1/4‘PRFD2 
2 • S  3  >  RAD 


S22-RAD2+1/4*PRFD 

s2*prpd 


>] 

)1 
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t,  ■  arccos 
4 


Sj2-RAD2+1/ 4 • PRFD4 

s3*prfd 


where  the  maxima  and  minima  are  taken  to  ensure  values  of  tr  and  0, 
respectively,  for  in  degenerate  triangles. 

By  definition  of  cosine  we  have 


with  the  minimum  taken  to  guarantee  a  value  of  zero  should  triangle 
RA2A3  be  degenerate. 

Should  the  burnout  criteria  fail,  the  sliver  is  not  yet  extinguished. 
Then  we  calculate 


area  triangle 

A^AjP  s> 

1/2 

•  RAD  • 

S2  * 

sin 

(xp 

area  triangle 

AiA2q  = 

1/2 

•  RAD  • 

S3  * 

sin 

Cx2) 

area  triangle 

= 

1/4 

•  PRFD* 

S1  * 

sin 

C°) 

2 

area  sector  A^PQ  =  1/2*  (AjAjAj-t^Tj)  *RAD 

2 

area  sector  A3PR  *  1/8*  (Tj-cr-A^A^ -PRFD 

area  sector  A2QR  =  1/8* (t^-o-AjA^A^) *PRFD2 
Then  the  sliver's  end  area,  PQR,  is  given  by 
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E  =  area  A^A^P  +  area  A^PQ  +  area  A^A2Q 

-  area  A^  Vj  “  area  A2AjR  -  area  AjPR  -  area  A2QR 

or,  more  simply, 

E  ■  1/2  RAD*  S2  •  sin(t^)  +  RAD*  (A3A^A2-Tj-T2)+S3*sin(T2) 

-  area  A^A2A3 

-  1/4*PRFD*  jsj-sinO)  +  1/2* PRFD*  (t j+t^o-A^Aj-AjA.^) 

and  its  contribution,  VOL,  to  grain  volume  is 
VOL  =  E-GRL 

Next,  note  that 

length  arc  PQ  ■  RAD  •  (A^Aj  -  -  t2) 

length  arc  PR  =  1/2  •  PRFD  •  (j3  -  o  -A^A^ 

length  arc  QP  =  1/2  •  PRFD  •  (t4  -  a  -  AjA^) 

Then  the  sliver's  contribution,  SURF,  to  grain  surface  area  is  twice  its 
end  area,  E,  plus  its  lateral  surface  area,  or 

SURF  =  2  •  E  + 

GRL  •  j  RAD  •  (A2AjA3  -  -  t2)  +  1/2  PRFD  •  (tj  +  t4  -  2a  -AjA^-A^A^ 
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When  extinction  occurs,  however,  we  have 


VOL  =  SURF  *  0 

C.  INNER  SLIVERS 


As  in  the  case  of  outer  slivers,  we  require  that  every  inner  sliver 
remain  within  its  original  boundaries  throughout  burning  until  extinction. 
Initially,  an  inner  sliver  appears  as  A^Aj  in  Figure  4  with  A.,  A2,  and 
A3  the  centers  of  the  sliver's  three  perforations. 

We  assume  lengths  S^,  82,  and  Si  to  be  known  and  calculate  angles 
Aj,  A2,  Aj  and  the  area  of  triangle  A.AjA^  in  the  same  way  as  in  the 
case  of  the  outer  sliver. 

As  burning  progresses,  the  sliver  erodes  as  shown  in  Figure  S  with, 
possibly,  some  degenerate  angles.  We  require  knowledge  of  the  grain's 
instantaneous  length,  GRL,  and  perforation  diameter,  PRFD.  As  before, 
the  surfaces  of  all  perforations  are  assumed  to  erode  at  the  same  rate. 
Burnout  of  the  sliver  occurs  when  either 
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ir/2 


or 


GRL  <  0 


An  error  condition  occurs  when  an  expanding  perforation  crosses  any 
portion  of  the  line-of-centers  c£  the  two  opposite  perforations  before 
either  of  them  engulfs  that  point,  signaling  violation  of  our  requirement 
that  the  sliver  remain  within  its  original  boundaries. 

The  critical  points  are  clearly  the  feet  of  altitudes  which  cut 
opposite  sides  internally,  as  in  Figure  6.  The  error  condition  arises 
whenever 


dj  <  min  (d2>  dj) 

This  is  equivalent  to  having  the  altitude  be  less  than  the  smaller  of 
the  two  segments  into  which  its  foot  divides  the  opposite  side.  In 
Figure  7,  suppose 
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h  <  min  (a,b) 


h  <  a,  h  <  b 


a  <  y  ,  fl  <  5 


a  <  ir/4,  0  <  ir/4 


and  conversely.  Thus,  the  error  arises  only  in  obtuse  triangles  where 
both  acute  angles  are  smaller  than  tt/4.  We  note  that  the  inner  sliver 
error  condition  need  only  be  checked  once,  based  as  it  is  on  original 
grain  dimensions  alone. 

Should  the  error  arise  due  to  a  problem  with  an  inner  sliver,  the 
grain  may  be  repartitioned  to  attempt  to  cure  the  fault.  A  good  policy 
is  to  try  to  avoid  grossly  obtuse  slivers. 

Referring  once  more  to  Figure  5,  we  have  by  definition  of  cosine 


■  arccos 


*  arccos 


s  arccos 


J^min  ^1,  PRPD)J 

[nin  (l-  pire)] 


f*ln  prfd)] 


where  the  minima  are  taken  to  ensure  a  value  of  zero  for  t . ,  in 
degenerate  triangles. 

Should  the  burnout  criteria  fail,  the  sliver  is  not  yet  extinguished. 
Then  wo  calculate 


area  triangle  AjAjR  =  1/4  •  PRFD  •  S^sin  (t19) 
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area  triangle  A^AjQ  ■  1/4  •  PRFD  •  S2  •  sin  (t13) 

area  triangle  AjAjP  *  1/4  ♦  PRFD  •  •  sin  23^ 

2 

area  sector  AjQR  ■  1/8  •  PRFD  •  (AjAjAj  -  -  Tj2) 

2 

area  sector  AjPR  ■  1/8  •  PRFD  *  (AjAjAj  “  tj2  "  T23^ 

2 

AJ'ci  sector  AjPQ  *  1/8  •  PRFD  •  (AjAjA2  *■  -  T23^ 

Then  the  sliver's  end  area,  PQR,  is  given  by 

E  ■  area  A^Aj  -  area  A^R  -  area  AjA^Q 

-  area  A^P  -  area  A^QR  -  area  A2PR  -  area  A3PQ 

or,  more  simply, 

E  ■  area  AjAjAj 

-  1/4  •  PRFD  •  |  S3*sin(t12)  +  S2*sin(T13)  +  S^sinCr^) 

+  PRFD  •  (it/ 2  -  x12  -  t13  -  t23) 

and  its  contribution,  VOL,  to  grain  volume  is 
VOL  ■  E*GRL 
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Next,  note  that 


length  arc  PQ  -  1/2  •  PRFD 
length  arc  QR  ■  1/2  *  PRFD 
length  arc  PR  *  1/2  •  PRFD 


•  -  t13  -  t2;}) 

•  (A2A1Aj  -  t12  -  t13) 

•  CAiA2A3  "  t12  “  t23‘) 


Then  the  sliver's  contribution,  SURF,  to  grain  surface  area  is  twice 
its  end  area,  E,  plus  its  lateral  surface  area,  or 

SURF  »  2  ■  E  +  GRL‘PRFD‘  (n/2  -  T12  -  -  t23) 


When  extinction  occurs,  however,  we  have 


VOL  «  SURF  -  0 


III.  SUBROUTINES  "GENIS"  AND  "GENOS" 

Subroutines  GENIS  (GENERAL  INNER  SLIVER)  and  GENOS  (GENERAL  OUTER 
SLIVER)  employing  the  above  calculational  methods  have  been  written  in 
FORTRAN  and  are  listed  in  the  Appendix. 

For  each  class  of  slivers,  two  arrays  are  necessary,  S(3)  and  A(4). 
In  array  S,  the  calling  routine  must  store  the  three  sides  of  triangle 
A1A2A3  (Figures  2  and  4)  taking  care,  in  the  case  of  the  outer  sliver, 
that  8(1)  contains  the  chord  segment  joining  the  centers  of  the  two 
perforations  bounding  that  sliver.  Array  A  is  a  work  array  which  will 
be  used  to  store  the  angles  of  AjA2A3  and  its  area  for  later  use. 

PRFD  is  the  grain's  instantaneous  perforation  diameter,  GRL  its 
instantaneous  length  and,  for  outer  slivers,  RAD,  the  grain's  radius. 
Units  of  measurement  should,  of  course,  be  uniform. 

The  routino  will  return  SURF,  the  sliver's  contribution  to  grain 
surface  area,  and  VOL,  its  volume  contribution. 

Branching  within  each  routine  is  controlled  by  argument  IFLAG 
which  the  calling  routine  should  initialize  to  a  negative  number  and  not 
touch  thereafter.  Extinction  of  the  sliver  will  be  signaled  by  a  value 
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of  IFLAG  «  1.  Thus  complete  extinction  of  the  grain  occurs  when  the 
branch  flag  for  each  class  of  slivers  is  1. 

IV.  A  SAMPLE  CASE 

Suppose  that  a  nineteen-perforated  grain,  one-sixth  of  which  is 
shown  in  Figure  8,  is  characterized  by  three  independent  webs  WI,WM,WO 
and  some  perforation  diameter,  PD,  and  length  GL. 

There  are  four  classes  of  slivers: 

(i)  6  inner  slivers  congruent  to  FED 

(ii)  6  inner  slivers  congruent  to  BED 

(iii)  12  inner  slivers  congruent  to  BEC 

(iv)  12  outer  slivers  congruent  to  CIHB. 

We  note  that  grain  diameter,  D,  is  given  by 

D  ■  5*PD  +  2  •  (WI  +  WM  +  WO) 
and  we  have  segments 

FE  -  WI  +  PD 

EC  ■  WM  +  PD 

FC-FE+EC-W1  +  WM  +  2-PD 

AB  “  1/2  •  FC"  <*  1/2  •  (WI  +  WM)  ♦  PD 

Fi  =■  AB  •  /T  »  {  1/2  *  (WI+WM)  +  PD}  •  /f 

EB  «^BJ  2  +  EJ  2  =  l[(¥b  -  F7)2  +  (1/2 -FE) 2 

-  1/2  •  1(f.(WM+PD)2  +  (WI+PD)2  , 
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Let  arrays  SIDE  (3,4)  and  ANGL  (4,4)  be  used  for  our  side  and  work  arrays, 
respectively.  Then,  considering  each  sliver  class  in  turn,  set: 

(i)  SIDE  (1,1)  =  SIDE  (2,1)  »  SIDE  (3,1)  ■  FE  »  WI  +  PD 

(ii)  SIDE  (1,2)  -  SIDE  (2,2)  -  Iff  -  1/2  •  ^3(WM+PD)2  +  (WI+PD)2 

SIDE  (3,2)  -  Iff  -  WI  +  PD 

(iii)  SIDE  (1,3)  ■  Iff  ■  1/2  •  (WI+WM)  •  PD 

SIDE  (2,3)  ■  Iff  ■  1/2  •  3 (WM+PD) 2  +  (WI+PD)2 

SIDE  (3,3)  -  EC  -  WM  +  PD 

(iv)  SIDE  (1,4)  -  BC  -  1/2  •  (WI  +  WM)  +  PD 

SIDE  (2,4)  -  PC  »  WI  +  WM  +  2  •  PD 

SIDE  (3,4)  -  FB  -  (1/2  •  (WI+WM)  +  PD)  ■  /T 

where,  as  noted,  SIDE  (1,4)  must  be  the  indicated  segment,  BC,  for  the 
outer  slivers. 

Let  array  NCHECK  (4)  be  set  aside  for  the  branch  flags. 

Now,  while  the  subroutines  can  be  called  throughout  the  burning 
of  the  grain,  it  is  efficient  to  employ  the  more  obvious  method  of 
subtracting  the  nineteen  cylindrical  perforations  from  the  grain  to 
deduce  surface  area  and  volume  until  at  least  one  sot  of  perforations 
touch.  That  is,  until  twice  the  depth  burned  is  WEB,  where 

WEB  =  min  (WO,  WM,  WI,  BC  -  PD,  BE  -  PD,  GL) 

=>  min  (WO,  WM,  WI,  1/2  •  (WI+WM), 

1/2  •  ^3« (WM+PD)2  +  (WI+PD)2  -  PD,  GL) 
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Afterward,  we  employ  GENIS  and  GENOS,  remembering  to  sum  the  individual 
sliver  contributions  to  produce  SSUM  and  VSUM,  the  grain's  surface  area 
and  volume,  respectively,  at  any  instant, 

A  rendering  of  a  program  which  determines  grain  surface  area  and 
volume  for  our  burning  grain  with  a  fixed  depth-burned  increment,  DELTA, 
is  in  the  Appendix.  It  solicits  grain  dimensions  from  the  user  and 
simply  prints  a  table  of  the  surface  and  volume  for  the  grain  at  increas¬ 
ing  depths  burned  through  extinction.  The  ACOS  function  is  included 
because  the  FORTRAN  library  of  the  minicomputer's  operating  system  doesn't 
contain  it.  The  outer  grain  surface  is  assumed  to  erode  at  the  same  rate 
as  the  perforation  surfaces  although  this  need  not  have  been  forced. 

Further  development  is  possible.  By  adding  a  third  subscript  to  our 
SIDE  and  ANGL  arrays,  the  program  could  be  extended  to  handle  several 
such  grains.  If  the  intermediate  results  were  saved,  it  would  be 
possible  to  present  them  in  plotted  form.  This  last  has  been  done  in 
Figures  9  and  10  which  display  the  ratio  of  surface  area  to  initial 
surface  area  plotted  against  fraction  of  grain  burned  for  two  hypothetical 
grains.  The  first  has  all  three  webs  equal,  while  the  second,  Figure  10, 
has  three  unequal  webs, 
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Figure  9.  Surface  Area  Plot  -  Grain  with  Equal  Webs 
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SAMPLE  19-PF  PROPELLANT  WITH  DIFFERENT  WEBS 


Figure  la  Surface  Area  Plot  -  Grain  with  Unequal  Webs 
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APPENDIX 

LISTING  OF  SAMPLE  PROGRAM  FOR  NINBTBBN- PERFORATED  GRAIN  PLUS  LISTINGS 
OF  SUBROUTINES  "GENOS"  AND  "GEN IS", 
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FORTRAN  IV 


V01C-03G  THU  13-FEB-79 


PAGE  001 


C 

SAMPLE  PROGRAM  EMPLOYING  #GENOS*  AND  *GEN!S* 

c 

TO  PRODUCE  SURFACE  t  VOLUME  FOR  19 -PERFORATED 

c 

p 

GRAIN  UITH  THREE  INDEPENDENT  WEBS. 

0001 

L 

DIMENSION  SIDE<3#4)*ANGL<4*4) *NCHECK<4) 

0002 

p 

DATA  NCHECK/4#-!/,  lDONE/0/*PI/3. 14139263/ 

0603 

L 

TYPE  1030 

0004 

ACCEPT  1000 *UI*UM*UQ* PD* GL* DELTA 

0000 

p 

D-5 .  xtPD+2 .  #(UI+UM+UO) 

0006 

L 

SIDEU.D-UI+PD 

01107 

SIDEC2* 1) -BIDE(1* 1) 

0008 

p 

SIDE<3* 1)  -SIDEU* 1) 

0009 

L 

SIDE(1.2>-.5*SQRT<3.*<UM+PD>**2  +  (UI+PD)#*2) 

0010 

SIDE(2»2) "SIDE( 1*2) 

0011 

p 

SIDG<3*2) “SIDE(1* 1) 

0012 

L 

SIDE(1*3)«.3#<UIH]MHPD 

0013 

SIDE<2*3) »SIDE< 1*2) 

0814 

p 

SIDG(3*3) -UM+PD 

0015 

SIDE(1*4) -SItEd.3) 

0016 

SIDE<2*4) «2.*SIDEd»3) 

0017 

r 

8IDE(3*4)  -SlDEd*3)#SQRT<3.) 

0018 

U 

p 

UE8-AMINl(U0,UM,UI*.5*aJI+UM)*SIDE(l*2)-PD*GL) 

0019 

L 

TYPE  1010*1)1*  DM*  DO*  PD*  6L*  DELTA 

0020 

JSTEP-0 

0021 

100 

DEPTH«JSTEP*DELTA 

0022 

PRFD«PD+2.*DEPTH 

0023 

GRL-GL-2.#DEPTH 

0024 

RAD«.S#D-DEPTH 

0025 

p 

IF<2.*DEPTH,GT.UEB)  GOTO  200 

0027 

L 

TA«PI#(RAD#*2-4.75>WPRFD**2) 

0028 

SSUM-2 .  *TA+GRL*P  I*(2 ,  #RAD+19.#PRFD) 

0029 

VSUM-TA*GRL 

0030 

p 

GOTO  300 

0031 

200 

I  DONE a0 

0032 

SSUM-0, 

0033 

p 

VSUM-0. 

0034 

L 

DO  250  1-1*2 

0035 

CALL  GEN  I S ( S I DE ( 1 , 1 > » ANGL  C 1  *  I ) *  PRPD . GRL *  NCHECK ( I ) * ' TA *  TV) 

0036 

SSUM-S8UM+6 . *TA 

0037 

VSUM-VSUM+6 . #TV 

0038 

1  DUNE  - I  DONE UNCHECK  ( I ) 
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FORTRAN  IV 


V01C-03G  THU  13-FEB-79 


PAGE  002 


0839  250  CONTINUE 

0040  CALL  GEN  IS  £S I DE  £  1 . 3) . ANGL (1.3).  PRFD.  GRL. NCHECK (3) .  TA.  TV) 

0041  SSUM«SSUM+t2.#TA 

0042  VSUM-V8UM+12 . #TV 

C 

0043  CALL  GENOS £8 IDE (  l , 4) * ANGL (U 4)# PRFD.GRL.RAD. NCHECK ( 4) . TA. TV) 

0044  SSUM«SSUM+12.#TA 

0045  VSUM-VSUM+12.#TV 

0046  IDONE- IDONE+NCHECK (3) +NCHECK  £4) 

C 

0047  300  TYPE  1020. DEPTH, SSUM.VSUM 
8048  IF C IDONE. EQ. 4)  GOTO  400 

0050  JSTEP-JSTEP+l 

0051  GOTO  100 

C 

0052  400  STOP  'GRAIN  EXTINGUISHED!' 

C 

0053  1000  FORMAT! 10E 12.0) 

0054  1010  FORMAT!/"  INNER  UEB  -  '.F10.4 
«  /'  MIDDLE  UEB-  '.FIB. 4 

*  /'  OUTER  UEB  •  '.FIB. 4 

«  /"  PERF  DIAM  -  '.FIB. 4 

*  /*  GRAIN  LEN  •  '.FIB. 4 

*  /'  BURN  I NCR  ■  '.FIB. 4 

$  //'  DEPTH  BURNED  SURFACE  AREA  VOLUME'/) 

0055  1020  FORMAT! IX. FI 0.4. 5X.F 10.4. 5X» FI 0.4) 

0056  1030  FORMAT!'  INNER. MIDDLE. OUTER  UEBS.PERF  DIAM'. 

$  '.GRAIN  LENGTH. DEPTH-BURNED  INCREMENT  ?') 

END 


0057 


FORTRAN  IV  V01C-03Q  THU  1S-FEB-F9 

0001  FUNCTION  ACOS(U) 

0002  AC08»ATAN2  <  SORT  ( i .  -U>W<2)  » U) 

0003  RETURN 

0004  END 
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FORTRAN  IV  V01C-03G  THU  15-FEB-79  PAGE  001 

0001  SUBROUTINE  GEN0S<3* A*PRFD*GRL* RAD. IFLAG. SURF. VOL) 

C 

C  SUBROUTINE  *GENOS*i  CALCULATE  SURFACE  AREA  AND  VOLUME  FOR  A 
C  GENERAL  OUTER  SLIVER  OF  A  BURNING  GRAIN 

C  WITH  LENGTH  «  GRL,  RADIUS  •  RAD*  AND 

C  PERF  DIAM  -  PRFD. 

C 

0002  DIMENSION  S(3).A<4> 

C 

0003  IF(IFLAG)  10*20*30 

C 

C  INITIAL  PASSi  IFLAG  HAS  SET  NEGATIVE  BY  CALLING  ROUTINE. 

C  STORE  ANGLES  A1*A2*A3  AND  AREA  OF  TRIANGLE 

C  WITH  SIDES  S(1)*8(2)*S(3)  INTO  AC1>* .. ,*A<4> 

C 

0004  10  All)  -  ACQS((S<2)^24€(5)**2-Sa>**2>/(2.*S(2)*S(3))) 

0005  A ( 2 )  ■  AC0Sasa)**2+SI3)**2-S(2)>K#2)/(2.n«sn)>KS<3))) 

0006  A(3)  •  ACOS(  (S<  1)##2+B(2)#>k2-S(3)#*2)/<2.#S(  1  )*S<2) ) ) 

C 

0007  A(4)  -  .5*SU)*S<3)*S!N<A(2>> 

C 

C  ...AND  SET  FLAG  TO  ZERO  TO  BYPASS  INITIALIZATION  HEREAFTER. 

C 

0008  IFLAG  -  0 

C 

C  SUCCEEDING  PASSES  UNTIL  BURNOUT!  FIRST  DETERMINE  AUXILIARY  ANGLES 
C 

0009  20  TAU1  ■  ACOS (AMIN  1 <  1 . *  <SC2)H<#2+RAD>!<#2-.25#PRFD»W<2)/'(2.#S(2)#RAD) )) 

0010  TAU2  ■  ACOS < AM IN  1 < 1 . * (S ( 3) •WCHRAD##2- . 25*nPRFDHot<2 ) / C 2 . *S ( 3) #RAD ) > ) 

O0U  TAU3  ■  AC0B(AMAXl(M.*<S<2)*#2H*AD##2+.25^RFD*#2:>/(S<2>#PRFDm 

0012  TAU4  -  AC0S(AMAX1 (-1 . * (S  (3)##2-RAD##2+.25*PRFD>M‘2)/(S(3)»#PRFD) ) ) 

C 

0013  SIG  -  ACOS (AMIN  1 ( 1 . *S(1)^PRFD) ) 

C 

C  ...THEN  CHECK  ERROR  CONDITIONS... 

C 

0014  IF(TAU3.LT.A(3)  .OR.  TAU4.LT.A(2))  STOP  '*GENOS#  ERROR!* 

C 

C  ...IF  OK*  BRANCH  IF  SLIVER  FAILS  BURNOUT  CRITERIA. 

C 

0016  IF(TAUl+TAU2.LT.A< 1)  .AND.  GRL.GT.0. )  GOTO  25 

C 

C  SLIVER  JUST  BURNED  OlJTi  SET  FLAG  TO  BYPASS  AREA  t  VOLUME  CALCULATIONS. 
C 

0018  IFLAG  -  1 

0019  GOTO  30 

C 

C  SLIVER  NOT  BURNED  OUT!  DETERMINE  END  AREA*  VOLUME,  AND  SURFACE  AREA 
C 

0020  25  E  •  .5HflAD#(S(2)#SIN<TAUl>+RAD*<A(l>-TAUl-TAU2>+S(3>#SIN<TAU2)) 

*  ~A(4)-.25#PRFD#<3(l)#SIN(SIG)+.5*PRFD#(TAU3+TAU4-2.#SIG 
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FORTRAN  IV  V01C-03G  THU  15-FEB-?9  PAGE  002 

*  -A<2)-A<3) ) ) 

C 

0021  VOL  «  E#GRL 

C 

0022  SURF  ■  2.^4GRL#<RAD#<A(1>-TAU1-TAU2)+.5*PRFD#CTAU3+TAU4“2.#SIG 

*  -A(2)-A<3) ) ) 

C 

C  ...AND  RETURN. 

C 

0023  RETURN 

C 

C  SLIVER  18  BURNED  OUTi  RETURN  WITH  ZERO  VOLUME  AND  SURFACE  AREA. 

C 

0024  30  VOL  -  0. 

002S  SURF  -  0. 

0026  RETURN 

C 

002?  END 


FORTRAN  IV 


V01C-03G  THU  15-FEB-79 


PAGE  00 L 


0001  SUBROUTINE  GENISCS.A.PRFD.GRL, IFLAG. SURF, VOL) 

C 

C  SUBROUTINE  *GENl8*i  CALCULATE  SURFACE  AREA  AND  VOLUME  FOR  A 
C  GENERAL  INNER  SLIVER  OF  A  BURNING  GRAIN 

C  WITH  LENGTH  ■  GRL  t  PERF  DIAM  -  PRFD. 

C 

0002  DIMENSION  8(3). A(4) 

0003  DATA  P12/  1 .9707963  / 

C 

0004  IF(IFLAG)  10.20.30 

C 

C  INITIAL  PASSt  IFLAG  UAS  SET  NEGATIVE  0Y  CALLING  ROUTINE. 

C  STORE  ANGLES  Al.A2.A3  AND  AREA  OF  TRIANGLE 

C  UITH  SIDES  S( 1) .8(2) *S(3)  INTO  A(l>. . . ,.A(4) 

C 

0005  10  Ad)  •  AC0S(<S<2)#*2+B<3)##2~Sd)##2)/'<2.#S<2)#S(3))) 

0006  ACS)  «  AC0S(CS(l)**2+SC3)##2-S<2)*#2)/(2.#S(t>*S<3))) 

0007  A(3)  -  AC0S((S(n*#2+S(2>##2“S<3)*#2)/<2.#S<l)*5(2))> 

C 

0000  A(4)  -  .5#S(1)#S(3)*SIN(A(2>) 

C 

C  ...CHECK  FOR  ERROR  CONDITIONi  FIND  IF  TRIANGLE  ACCEPTABLE... 

C 

0009  J  -  0 

0010  DO  15  I  -  1.3 

0011  IF(A(I),LT..3#P12)  J  -  J+l 

0013  15  CONTINUE 

0014  IF(J.GT.l)  STOP  '#GENIS#  ERROR  I* 

C 

C  IF  OK.  SET  FLAG  TO  ZERO  TO  BYPASS  INITIALIZATION  HEREAFTER. 

C 

0016  IFLAG  "  0 

C 

C  SUCCEEDING  PASSES  UNTIL  BURNOUT:  FIND  AUXILIARY  ANGLES 
C 

0017  20  TAU12  -  AC08(AMIN1U..S(3)/PRFD>) 

0010  TAU13  ■  AC0S(AMINl(l.*S(2)/'PRFD)) 

0019  TAU23  -  AC08 (AMIN  1 ( 1 . . S ( I ) /PRFD) ) 

C 

C  ...AND  BRANCH  IF  SLIVER  FAILS  BURNOUT  CRITERIA. 

C 

0020  I F < TAU 12+TAU 13+TAU23 , LT. P 12  .AND.  GRL.GT.0. )  GOTO  25 

C 

C  SLIVER  JUST  BURNED  OUT:  SET  FLAG  TO  BYPASS  AREA  «.  VOLUME  CALCULATIONS. 
C 

0022  IFLAG  -  1 

0023  GOTO  30 

C 

C  SLIVER  NOT  BURNED  OUT:  DETERMINE  END  AREA.  VOLUME.  AND  SURFACE  AREA 
C 

0024  25  E  -  A (4)  -  .25#PRFD»t<(S(  1  )WSIN(TAU23)+S<2)#8IN(TAU13) 
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*  4S<3)*8IN<TAU12)+FRFD#CP12-TAU12“TAU13-TAU23;> ) 

C 

0022  VOL  -  E*GRL 

C 

0026  SURF  ■  2 . *£+GRL#PRFD#(P I2-TAU12-TAU 13-TAU23) 

C 

C  ...AND  RETURN. 

C 

002?  RETURN 

C 

C  SLIVER  IS  BURNED  OUT>  RETURN  UITH  ZERO  VOLUME  AND  8URFACE  AREA. 
C 

0026  30  VOL  -  0. 

0029  SURF  -  0. 

0030  RETURN 

C 

0031  END 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below ;  teai  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1.  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/ contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6,  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 
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